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Abstract 

If axions or other low-mass pseudoscalars couple to electrons ("fine 
4^ structure constant" Ua) they are emitted from red giant stars by the 

Compton process 7 + e ^ e + a and by bremsstrahlung e + (Z, A) 
iZ,A) + e + a. We construct a simple analytic expression for the 
energy-loss rate for all conditions relevant for a red giant and include 
axion losses in evolutionary calculations from the main sequence to 
the helium flash. We find that Ua ^ 0.5 xlO"^*^ or ma ^ 9 meV/ cos^ f3 
lest the red giant core at helium ignition exceed its standard mass 
by more than 0.025 A^©, in conflict with observational evidence. Our 
bound is the most restrictive limit on cIq, but it does not exclude the 
possibility that axion emission contributes significantly to the cooling 
of ZZ Ceti stars such as G117-B15A for which the period decrease 
was recently measured. 
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1 Introduction 



The cooling rate of the ZZ Ceti star G117-B15A as determined from the 
decrease of its pulsation period appears to be somewhat faster than can be 
accounted for by standard photon cooling. Isern, Hernandez and Garcia- 
Berro [1] speculated that this discrepancy was an indication for a novel cool- 
ing agent, notably for the emission of "invisible axions". 

Axions [2] are low-mass pseudoscalar particles that couple to electrons by 
virtue of a Lagrangian density 

Ant = -igi^el^i^ed (1) 

where is a dimensionless coupling constant, r/jg is the electron Dirac field, 
and a the axion field. We shall also use the "axionic fine structure constant" 

eta = 5'^/47r and = cta/10~^^. (2) 

In a certain class of models (DFSZ axions) the Yukawa coupling is 

5 = 2.8 X lO"^"* m^ev cos^ /3 (3) 

where cos^ /3 is a model-dependent parameter which we shall always set equal 
to unity, and mmeV is the axion mass in units of 1 meV = 10""^ eV. Then, 

a26 = 0.64x10-2 m^^v- 

The main energy-loss mechanism in a white dwarf is bremsstrahlung emis- 
sion e+(Z, A) — )• (Z, A) + e+a. Isern, Hernandez and Garcia-Berro [1] favored 
an axion mass of 8.4 meV, equivalent to = 0.45, in order to explain the 
cooling rate of G117-B15A. 

Of course, this interpretation is very speculative and so, naturally one 
wants to know if it is consistent with other astrophysical phenomena that 
might be affected by axion emission. For example, the overall white dwarf 
luminosity function leads to a constraint of ~ 1-0 [3]. 

Another constraint was derived by Wang [4] who required that axion 
cooling would not prevent carbon ignition in accreting white dwarfs so that 
type I supernova explosions can occur. Wang's bound, based on a simple 
analytic estimate, is a26 ~ 6 or ^ 30meV. 

Horizontal-branch stars have a nondegenerate, helium-burning core which 
would emit axions dominantly by the Compton process 7 + e— )-e + a. A 
crude bound is based on the requirement that the energy-loss rate should 
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not exceed 100 erg g~^ or else the HB lifetime would be reduced by more 
than about a factor of two, in conflict with the observed number of HB stars 
in globular clusters. Then one finds the bound ~ 5 [5]. 

The potentially most restrictive argument discussed in the literature was 
put forth by Dearborn, Schramm and Steigman [6]. They considered the 
impact of axion emission on red giants near the helium flash; for ~ 0.16 
they found helium ignition to be suppressed entirely which would clearly 
contradict the mere existence of the horizontal and asymptotic giant branches 
observed in stellar systems. Unfortunately, they used emission rates which 
did not take degeneracy effects properly into account; near the center of a 
red giant before helium ignition they overestimate the energy-loss rate by 
as much as a factor of 10 (see below). Still, their adjusted limit on a^^ is 
only a factor of 2 or 3 above the value favored to explain the cooling rate 
of G117-B15A, and so, it seems worthwhile to revisit the helium ignition 
argument with a more appropriate energy-loss rate. 

2 Energy-Loss Rate 
2.1 Compton Process 

The simplest possibility to produce axions by virtue of their coupling to 
electrons is the Compton process 7 + e — )• e + a [7]. In the nonrelativistic 
limit one finds a cross section a = Airaaaio^ j'im'l with a = 1/137 and u the 
photon energy. A simple integral over the initial-state photon phase space 
then yields the energy-loss rate per unit mass 



where ~ 1.0173, is the number of electrons per baryon, mjv is the nu- 
cleon mass which is used for an approximate conversion between the number 
density of baryons and the mass density of the medium, and Tg, = T/IO* K. 

The factor F accounts for relativistic corrections as well as for degeneracy 
effects and the nontrivial photon dispersion in a medium. For our purposes, 
the most severe deviation from F = 1 occurs at the center of a red giant before 
the helium flash. Taking p = 10^ g/cm"^ and T = 10® K = 8.6 keV as nominal 
values, the plasma frequency is 18keV and the electron Fermi momentum 
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is 409 keV whence the degeneracy parameter is rj = (fi — me)/T = 16.7. 
Typical blackbody photons have an energy of about 3T whence corrections 
from a "photon mass" remain moderate. Also, relativistic corrections to 
the emission rate are only about a 30% effect (Fukugita, Watamura and 
Yoshimura [8]). 

These authors also gave a table for F on a grid of T and p. For a fixed 
temperature, their values for F slightly increase with increasing density, con- 
trary to the expectation that degeneracy effects should decrease the emission 
rate. Upon closer scrutiny we are unable to find a Pauli- bio eking factor in 
their expressions of the phase-space integrals. We believe that the Comp- 
ton process must be suppressed by electron degeneracy which implies that 
bremsstrahlung dominates (see below). Therefore, a precise calculation for 
the degenerate regime is not warranted. In order to interpolate between 
degenerate and nondegenerate conditions, however, an estimate of the sup- 
pression factor Fdeg is useful. 

In the nonrelativistic limit electron recoils can be neglected so that the 
initial- and final-states have the same momentum. Therefore, F^eg is the 
Pauli blocking factor, averaged over all electrons, 

_ 1 r 2£p 1 / 1\ 

where p is the electron chemical potential and Ue the electron density. Then, 
1 Z""" „ ,„ 



Fdeg = pEdE yKj) 

where x = [E — fi)/T. For degenerate conditions the integrand is strongly 
peaked near x = (the Fermi surface) so that one may replace p and E with 
Pp and E'p, respectively, and one may extend the lower limit of integration 
to — oo. The integral then yields T so that 

Fdeg = 3EfT/pI, (7) 

where Ue = Pp/^^^ used. A Fermi momentum pp = 409 keV implies 
Ep = 655 keV; with T = 8.6 keV this gives F^eg = 0.10. Of course, there are 
relativistic corrections to this result. 
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2.2 Nondegenerate Bremsstrahlung 



The nondegenerate bremsstrahlung rate e + (Z, A) — )• (Z, A) + e + a was first 
calculated by Krauss, Moody and Wilczek [9] and e + e — )• e + e + a was 
added by Raffelt [10]. Ignoring screening effects which are a small correction 
for nondegenerate conditions, and allowing a chemical composition of only 
hydrogen (mass fraction X) and helium (mass fraction 1 — X) the energy-loss 
rate is 
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where Ts = T/IO* K as before and = g cm 



2.3 Degenerate Bremsstrahlung 

The degenerate bremsstrahlung rate was calculated in order to derive a bound 
on eta from white dwarf cooling times [3]. In this case screening effects must 
be included; otherwise the emission rate diverges. As a screening scale the 
electron Thomas-Fermi wave number was used, a common but incorrect prac- 
tice, which leads to an underestimate of the screening suppression because 
the main contribution is from the ions. Of course, because the screening scale 
enters logarithmically the error remains moderate — a factor of 2 or 3 for the 
white dwarf cooling rate. 

The axion emission rate for very degenerate matter relevant for white 
dwarfs and the crust of neutron stars was also calculated [11]. The main 
point was to include ion correlations in a strongly coupled plasma, a condition 
quantified by the parameter 

0.2275 ' ' 



aT Ts \A 

where Z is the charge of the ions, A their atomic mass, and n their density 
which determines the ion-sphere radius a = (3/47rn)^/"^. For F > 178 the 
ions arrange themselves in a bcc lattice while for F ^ 1 their correlations are 
weak. In a white dwarf F is typically between 20 and 150. 
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However, red giants near helium ignition are hot; for our standard set 
of parameters we find F = 0.57 which implies that Debye screening is still 
a reasonable description of the ion correlations. The electrons contribute 
little to screening because the Thomas-Fermi wave number is much smaller 
than the Debye scale; otherwise the plasma would not be degenerate. (For 
our standard red giant conditions the Thomas-Fermi wave number is about 
50keV while the Debye scale for the ions is 222 keV.) 

With these approximations one finds for the energy-loss rate [12] 
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with /3f = Pf / Ep the velocity at the Fermi surface. With ki) the Debye 
screening scale of the ions (density n, charge Ze) 
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For helium this is k,^ 



1 /3 

0.147 Pq /Tg. For our benchmark conditions we have 



/3| = 0.39 and then F = 1.7. 



2.4 Interpolation Formula 

The main region of interest to us is the degenerate red giant core. However, 
the hydrogen burning shell is entirely nondegenerate and also at a temper- 
ature of about 10® K so that a consistent treatment requires to implement 
axion emission everywhere in the star. To this end we interpolate between 
the degenerate and nondegenerate bremsstrahlung rates by 

e=(%D + e5')-'- (13) 

The nondegenerate Compton rate is switched off in the degenerate regime by 
means of a factor (1 + Fd"e|)"^^^ where Fdeg was given in Eq. (7). In Fig. 1 
we show the different rates as well as our interpolation as a function of p 
for T = 10® K. Interestingly, the total rate is nearly independent of density; 
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Figure 1: Axionic energy-loss rates for the processes discussed in the text for 
= 1,7" = 10^ K, and a composition of pure hehum. The soHd Hne is our 
interpolated emission rate. 

this is a coincidence at the given temperature because the Compton and 
degenerate bremsstrahlung rates vary with different powers of T. 

Dearborn, Schramm and Steigman [6] gave a table of their energy-loss 
rates. For a coupling constant = 1 and T = 10® K they used 20, 50, and 
201ergg~^s~^ at densities 10^, 10"*, and lO^g/cm"^. At the highest relevant 
density this is about a factor of 10 above our rate. 



In order to test the impact of axion emission on the evolution of red gi- 
ants we have included the interpolation formula described in the previous 
section in our stellar evolution code in analogy to our previous study of non- 
standard neutrino losses [13]. We have then calculated several evolutionary 
sequences from the main sequence to helium ignition with different axion 
coupling strengths a26- We used a chemical composition corresponding to 
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Mixture I of Ref. [13], i.e., to Z = and Yq = 0.239. The opacities were 

chosen for a Ross-AUer mixture; the impact of axion emission on the core 
mass is found to be the same for older Los Alamos ("AOL") [14] as well as the 
latest Livermore ("OPAL") [15] opacities, which have greatly improved the 
agreement between observations and stellar evolution theory in general. The 
mixing length parameter is taken to be 1.55. The plasma neutrino energy- 
loss rate was taken from Ref. [16]. The total stellar mass was 0.8 A^©; mass 
loss on the red giant branch was ignored. For other aspects of our stellar 
evolution calculations see Ref. [13] and references therein. 

We find that helium ignites at a core mass A^ig which is increased by 
the a26-dependent amount which is given in Tab. 1 and shown in Fig. 2. 
The coupling strength = 2 corresponds approximately to the case where 
helium ignition was suppressed in the calculations of Dearborn, Schramm 
and Steigman [6] if one corrects for the overestimate of their emission rate. 
Even for stronger couplings helium still ignites in our calculations, although 
for our largest value (a26 = 8) the core-mass increase is so large that, had 
we included mass loss, the entire envelope could have been consumed before 
helium had a chance to ignite. 

Even though our calculations do not reproduce the suppression of helium 
ignition, which is an overly conservative criterion to constrain axion emission, 
we believe that the core mass increase alone can be used to derive a significant 
limit on a26- 

Table 1: Increase of the core mass at helium ignition. 



"26 -^A^ig [Mq] 

0.0 0.000 

0.5 0.022 

1.0 0.036 

2.0 0.056 

4.0 0.080 

8.0 0.111 
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Figure 2: Increase of the core mass of a red giant at helium ignition due to 
axion emission. 

4 Discussion and Summary 

It was previously shown [17] that observations of globular cluster stars and of 
field RR Lyrae stars confirm the standard core mass at helium ignition A^ig 
to within about 5%, i.e., to within about 0.025 A^©. The main observational 
constraint is the maximum brightness reached by red giants, and the observed 
brightness of field RR Lyrae stars. We have previously used this method to 
constrain neutrino magnetic dipole moments [13]. 

A^ig depends slightly on the total stellar mass and on the chemical com- 
position (see [13, 17] for approximate analytic formulae); it is about O.49OM0 
for a helium content of 0.26 and a metallicity of 0.001. The systematic uncer- 
tainties of A^ig due to possible deviations of the opacities from a Ross-AUer 
metallicity mixture, due to the standard mixing length treatment of con- 
vection, mass loss on the red giant branch, and the numerical shell- shifting 
technique all seem to be much smaller than this limit [13, 18, 19]. 

A core-mass increase of 0.025 A^© corresponds approximately to a^^ = 
0.5, i.e., we find that globular cluster stars require that 

a„< 0.5x10-2^ or m„ < 9 meV/ cos^ /3. (14) 

This is the strongest bound currently available on the axion-electron coupling. 
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but it is not in conflict with the interpretation that axions could contribute 
significantly to the cooling of ZZ Ceti stars such as G117-B15A. 
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